Phase Boundary of the Boson Mott Insulator in a Rotating Optical Lattice 
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We consider the Bose-Hubbard model in a two dimensional rotating optical lattice and investigate 
the consequences of the effective magnetic field created by rotation. Using a Gutzwiller type vari- 
ational wavefunction, we find an analytical expression for the Mott insulator(MI)-Superfluid(SF) 
transition boundary in terms of the maximum eigenvalue of the Hofstadter butterfly. The de- 
pendence of phase boundary on the effective magnetic field is complex, reflecting the self-similar 
properties of the single particle energy spectrum. Finally, we argue that fractional quantum Hall 
phases exist close to the MI-SF transition boundaries, including MI states with particle densities 
greater than one. 
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Experiments on ultracold atoms in optical lattices 
opened up a new avenue to study correlated quantum 
states [![. The versatility of cold atom experiments hold 
promise for the experimental realization of many models 
that were first introduced for solid-state systems. 

One such model is the study of particles moving in 
a tight binding lattice under a magnetic field. When 
the magnetic flux per plaquette of the lattice becomes 
of the order of a flux quantum hc/e, the single particle 
energy spectrum forms a complicated self-similar struc- 
ture, known as the Hofstadter butterfly (Fig. [l} [2|]. It 
has not been possible to reach this regime in ordinary 
condensed matter experiments due to the required high 
magnetic fields. However, the ultracold atom experi- 
ments are extremely flexible and it should be possible 
to create required effective magnetic fields in optical lat- 
tice experiments. A conceptually simple way of creating 
an effective magnetic field is to rotate the optical lattice, 
as demonstrated in a recent experiment Q . Other means 
of creating effective magnetic fields have been discussed 
by a number of authors 0, [E Hi 0> HI • Although the re- 
cent demonstration of a rotating optical lattice was done 
for a shallow lattice, it should be possible to drive the 
system into the Mott insulator (MI) state by increasing 
the lattice depth. 

In this Letter, we study the Bose-Hubbard model un- 
der a magnetic field. Particularly, we consider a two di- 
mensional square lattice of spacing a with only nearest 
neighbor hopping. The magnetic field (or the effective 
magnetic field) strength can be expressed in terms of 
the dimensionless quantity <j>, which is the magnetic flux 
quantum per plaquette in the lattice (a 2 H/(hc/e), H be- 
ing the effective magnetic field). When the Landau gauge 
A = (0, Hx, 0) is chosen, the Hamiltonian for this system 
can be written as 

h = -tj24 a i eiAii + -^Y^fiiini-i) - (ij2™i, (!) 
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FIG. 1: (Color online) Maximum energy of the Hofstadter 
butterfly f(<f>) for a given <fr = p/q. This value is calculated as 
the maximum eigenvalue of the matrix A q — A.q(k x — 0, k y — 
0) (Eq. ©). 

where en (etj) is the bosonic annihilation (creation) op- 
erator at site i and hi = ajoj is the number operator. 
The tunnelling strength between nearest neighbor sites 
is given as t; U is the on-site interaction strength, and 
/i is the chemical potential. Magnetic field affects the 
Hamiltonian through Ay which is equal to ±27rm0 , if i 
and j have the same x coordinate ma and is otherwise, 
while the sign is determined by the hopping direction. 

We first review some of the properties of the single par- 
ticle spectrum by setting/7 = 0. This problem was first 
discussed by Hofstadter 0. The energy spectrum is ob- 
tained through the following difference equation (known 
as Harper's equation): 

E 

Cm+i + Cm-i + 2cos(27rm0 - k y )c m = —c m , 

where c m are the expansion coefficients of the wavefunc- 
tion, which has plane wave behavior along y in accor- 
dance with the translational symmetry in this direction. 
If <p is a rational number p/q, the wavefunction satisfies 
the Bloch condition c m+q = exp(ik x q)c m as a result of 
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the symmetry under g-site translation in the x direction. 
The allowed energies are then found as the eigenvalues of 
the q x q tridiagonal matrix: 



-ik x q 



1 

1 2 cos(27rm</> — k y ) 1 
1 



V 



•(2) 



We call the matrix formed by setting k x = k y = in @ 
A q . The maximum eigenvalue of A q yields the maximum 
energy of the system for a given </>. We define this energy 
as /(</>), which is a continuous but not differentiable func- 
tion (Fig. [T). To prove that the maximum eigenvalue is 
obtained from A g , we investigate the characteristic equa- 
tion for the matrix @, which is of the following form: 



t 
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2cos(k x q) — 2cos(k y q) — 0. (3) 



Two pairs of (k x ,k y ), namely (0,0) and Jjt/q, it/q) are 
sufficient to determine the band edges The (0,0) 
pair gives a smaller value for the k x and k y dependent 
terms. Since the E dependent part of ([3]) increases mono- 
tonically after a sufficiently large E, the greatest root is 
always obtained from the (0, 0) pair. 

We now turn to the interacting case with the dimen- 
sionless Hamiltonian: 

where t = t/U and // = fi/U are the scaled hopping 
strength and chemical potential. 

When the hopping term is dominant t > 1, one expects 
the system to be in a SF state, while in the opposite limit 
of strong interactions t <C 1, the system should go into 
the MI state. In this Letter, we investigate the transi- 
tion boundary between these two phases, and how this 



boundary is affected by the external magnetic field. The 
effect of the magnetic field on the transition boundary 
has been previously explored by strong coupling expan- 
sion for small magnetic fields by Niemeyer et. al. [lCl ]. 
and numerically within mean- field theory by Oktel et. al. 
[Tlj . Here we use a variational approach to provide an 
analytical expression for the transition boundary. 

We use a site dependent Gutzwiller ansatz to describe 
the system [12j. For the Bose- Hubbard model without 



magnetic field, this ansatz (and equivalent mean-field 
theory [H,[3|) gives an accurate description of the phase 
diagram. We introduce the variational wave function at 
each site I, 



\G)i = Ai\n a - 1); + |no)j + A[\n + 1) ; . 



(5) 



Since we investigate the behavior in the vicinity of the 
transition region, we consider small variations around the 
perfect MI state with exactly no particles per site, allow- 
ing for only one less or one more particle in a site. The 
variational parameters A; and A[ are assumed to be real, 
as complex A values can only increase the energy of the 
variational state. Total wavefunction is the direct prod- 
uct of these site wavefunctions |\&) = FJ^ |G};. Within 
the selected gauge, the magnetic Hamiltonian has trans- 
lational invariance in the y direction. The translational 
invariance in the x direction is broken by the magnetic 
field, but can be restored to a certain degree if the flux 
per plaquette is a rational number. Thus, taking (f> = p/q 
where p and q are relatively prime integers, the Hamil- 
tonian is invariant under translation by q sites in the x 
direction. This periodicity simplifies the calculation of 
the expectation value of the energy when we work with 
a supercell of 1 X q sites. Total wavefunction for such a 
supercell is \ty) s = YlIZo The expected value of the 
energy can then be written as follows 



WW) 



N s e, 



(0) 



where N s is the number of supercells. 

Keeping terms up to second order in the variational 
parameters A, the energy of a supercell is calculated as 
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- 2i{noA z A; +1 + ^n Q (n + 1)A Z AJ +1 + v^oK + i)A m Aj + (n + 1)A{A{ +1 
+ cos(2tt-Z) [n Af + 2^o(«o + 1)A ; AJ + (n + 1)(A|) 2 ] } 
+i [2(1 - n )Af + 27i o (A0 2 + n (n - 1)] + ^[A 2 - (A[) 2 - n ] 



(7) 



If the system favors to be in the Mott insulator state, the total energy of the system should be a minimum 
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FIG. 2: (Color online) The boundary of the Mott insulator 
phase for the first three Mott lobes. The figure is periodic in (f>. 
Magnetic field increases the critical value for t/U, as expected, 
however this increase is not monotonic. Transition boundary 
for two different values of fi/U are marked to display the 
complex structure of the surface. 



where all the variational parameters vanish. Thus, we 
can find the phase boundary as the point where the total 
energy ceases to be a local minimum in A. As a result, we 
demand that the matrix composed of the second deriva- 
tives of e with respect to the parameters (8 2 e/dAidAj, 
d 2 e/dA l dA' j , d 2 e/dA' i dA' J ) be positive definite, i.e. all 
eigenvalues be positive. This matrix can be written com- 
pactly as: 



F = -2t 



n A„ 



\fn (riQ + l)A q 

yjn (n + l)A g (n + l)A g 
2(1 - n + jl)l q 







2(n - p)l q 



where I q is the q x q identity matrix, and A 9 was intro- 
duced before (Eq. [2]). 

If we denote the eigenvalues and eigenvectors of A g by 
Ay and and those of F by A u and it, all X u can be 
expressed in terms of X v by taking 



bit 



due to the special block form of F. Then X u are obtained 
as: 



AJ = 1 - (1 + 2n )iX v =F yj [(1 + 2n a )tX u - l] 2 - 4{(n - /2)[1 - (n - fl)} -t(l + p)X v }. 

I 



The positive definiteness of F leads us to take A~ and set 
it to in order to determine the critical t value above 
which the perfect insulator state is destroyed. We find 
the boundary of the n^ 1 Mott lobe to be: 



t c = 



(no - p.)[l- (no - A)] 

(1 +£)/(#) 



(8) 



where uq — 1 < jl < hq. This boundary is plotted in 
Fig. for the first three Mott lobes. At <f) = 0, this 
formula reproduces the critical t value found in (Tj, [lj[ • 
Increasing magnetic field increases the critical hopping 
strength t c , however this increase is not monotonic. The 
complicated structure of the single particle problem is 
reflected in the transition boundary. Equation ([8]) is in 
excellent agreement with the numerical mean-field work 
& 

We can comment on the accuracy of our variational 
approach. Our result is exact within mean-field theory. 
At zero magnetic field the mean-field result for the tran- 
sition boundary is close to accurate Monte Carlo calcula- 
tions 15j, but it is not guaranteed that the mean-field 
description of the system would be valid under mag- 
netic field. Our variational wavefunction (and mean-field 
theory) disregards the correlations between fluctuations 



above the insulating state. Such correlations would be 
expected to wash out the fine structure of the transition 
boundary (Fig. 2). Nevertheless, one can expect a num- 
ber of features of the mean-field boundary to survive for 
the real system. The linear increase of the transition 
point for small magnetic fields, periodicity of the system 
with (j), and the central dip near = 0.5 should be qual- 
itatively correct. 

There is however one important way that the fluctua- 
tions around the Mott insulating state can become corre- 
lated. The Hamiltonian (|3|) supports bosonic fractional 
quantum Hall (FQH) states as discussed in a number of 
recent papers d, [H, [T3|. So far, such FQH states have 
been assumed to appear only in the region of low density 
where the number of particles per site is less then one. 
Here, we argue that states similar to bosonic FQH states 
should be present near the MI boundaries, even at higher 
densities. 

It is instructive to think about the behavior of the 
Hamiltonian for constant particle density by disregard- 
ing the last term. Let us assume that the particle density 
is equal to n = no + e, where no is an integer and e « 1 
is the decimal part of the density. With such incom- 
mensurate particle number, the system never goes into 
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FIG. 3: Schematic phase diagram near the ri^ Mott lobe. 
Dotted lines show the chemical potential as a function of hop- 
ping strength for systems with constant density (n) = no and 
(n) = no + e. FQH phases of "excess" particles, or holes are 
shown as the shaded regions. 

the MI state, but will always have a superfluid density. 
The chemical potential for this state, plotted on the jl, t 
plane, traces the outline of the Mott lobe as the interac- 
tion is increased (Fig. 3). However, if we think of the 
same system under a magnetic field that is commensu- 
rate with the excess particle density, another possibility 
presents itself. Specifically, considering a magnetic field 
so that <fi — 2e, it is possible for hq particles to form a MI 
state that is coexisting with a v = 1/2 Bosonic Laughlin 
state of the remaining e particles. At high enough inter- 
action, such a state would be preferable to a superfluid 
state as it avoids any interaction between the "excess" 
particles. The wavefunction of such a state can be ob- 
tained by symmetrizing the product of the Mott insulator 
state for no bosons with the v = 1/2 Bosonic Laughlin 
state for e particles. In general, separating the many par- 
ticle wavefunction into two parts and arguing that the 
overall properties can be deduced by thinking about the 
individual parts is not correct, as symmetrization may 
change the character of both parts considerably. In this 
case, however we can safely regard the excess particles as 
forming a correlated state above the Mott insulator, due 
to the full translational invariance of the MI state. One 
can write down an effective Hamiltonian for the excess 
particles. To the zeroth order, the change in the effective 
Hamiltonian would be just to replace t by (no + l)t, due 



to bosonic enhancement of the hopping. There will be 
higher order corrections to t and new non-contact inter- 
action terms between the excess particles due to fluctua- 
tions in the MI state. Such terms will be of higher order 
in (t/U), and can be neglected in the strongly interacting 
limit. One can also argue that as both the MI state and 
the Bosonic Laughlin state are gapped, it would not be 
energetically favorable to exchange particles between the 
two parts of the wavefunction. Similarly, one can argue 
that the overall state would be gapped in the strongly 
interacting limit. 

Treating such a state as a variational state, the energy 
difference from the MI state can be written as 

AE = {Uno -fx- t(n + l)f{4>)) e. (9) 

To first order in t/U , the term in parenthesis is the en- 
ergy needed to put one extra particle on to the Mott 
insulator. Thus, when it is favorable to put one extra 
particle on to the Mott state, it would be favorable to 
put more particles (up to e per site) and organize them 
into a FQH state. One can then expect the correlated 
state to exist within a band above the MI lobe (see Fig. 
[3]) . The same argument can be advanced for holes in a MI 
state, creating a FQH of holes below the Mott insulator. 
Experimentally these states would have distinct signa- 
tures appearing as extra steps near the MI steps of the 
Ziggurat structure of the trapped MI. Detailed proper- 
ties of these correlated states, as well as other correlated 
states near the transition boundary will be investigated 
elsewhere [l8j . 

In conclusion, we studied the phase boundary of the 
MI state of bosons in a rotating optical lattice. Using 
a Gutzwiller ansatz, we gave an analytical expression 
for the phase boundary in terms of the maximum en- 
ergy of the Hofstadter butterfly. We finally argued that 
analogues of FQH states will be found close to the MI- 
SF transition boundary including MI states with particle 
densities greater than one. 
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